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Abstract 

Early type massive stars drive thin, dense shells whose edges often show 
evidence of star-formation. The possibility of fragmentation of these shells, 
leading to the formation of putative star-forming clumps is examined with 
the aid of semi-analytic arguments. We also derive a mass-spectrum for 
clumps condensing out of these shells by performing Monte-Carlo simulations 
of the problem. By extending on results from our previous work on the 
stability of thin, dense shells, we argue that clump-mass estimated by other 
authors in the past, under a set of simplifying assumptions, are several orders 
of magnitude smaller than those calculated here. Using the expression for 
the fastest growing unstable mode in a shock-confined shell, we show that 
fragmentation of a typical shell can produce clumps with a typical mass > 10 3 
M . It is likely that such clumps could spawn a second generation of massive 
and/or intermediate-mass stars which could in turn, trigger the next cycle of 
star-formation. We suggest that the ratio of shell thickness-to-radius evolves 
only weakly with time. Calculations have been performed for stars of seven 
spectral types, ranging from Bl to 05. We separately consider the stability 
of supernova remnants. 
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1. Introduction 

Partial or full ring-like structures, which apparently are projections of 
shells, are often found in the interstellar medium (ISM). These dense shells 
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could possibly be driven by one of the numerous sources like, ionising ra- 
diation from young star-clusters, early type massive stars, blast-waves from 
supernovae, or energetic stellar winds. Detailed observations of these shells 
in various bands of the infrared wavelength have also revealed isolated sites 
of massive star formation (e.g. Deharveng, Zavagno & Caplan 2005). In 
the recent past a number of such sites have been reported, for instance in 
the HII region RCW79 (Zavagno et al. 2006), and RCW120 (Anderson et 
al. 2010). A catalogue of 600 such shells in the galactic disk was drawn up 
by Churchwell et al. (2006) as part of the GLIMPSE survey. The survey 
showed that a large proportion (~90 %) of these shells are thin, i.e. the shell 
thickness is less than a third of the outer shell radius, and driven primarily 
by massive stars (~ 86% shells), see also Zavagno et al. 2010. 

It is well-known that relatively high-mass stars emit powerful radiation 
that ionises gas in the local neighbourhood, and heats it to temperatures 
typically of the order of 10 4 Kelvin. This hot plasma propagates in the in- 
terstellar medium (ISM) at a highly supersonic speed (typical velocity of ex- 
pansion in the initial phases is ~ 10 3 km/s), whence it gradually equilibrates 
to a significantly lower temperature, of order a few thousand Kelvin. The 
expanding, roughly spherical volume of hot plasma, the so called Stroemgren 
sphere, cools primarily via collisional excitation of heavier elements such as 
Carbon, nitrogen, and oxygen, and sputtering of dust-grains while sweeping 
up a dense shell of gas in the ISM; the familiar snow-plough phase. This shell 
is confined by two shocks, first, due to the wind driving it, and second, due 
to the reverse shock resulting from the propagation of the shell in the ISM. 
The stability of such shells has been discussed by numerous authors, e.g. 
Elmegreen & Elmegreen (1978), Larson (1985), Elmegreen (1989), Vishniac 
(1983, 1994), Whitworth (1994), Wuensch & Palous (2001) and Anathpindika 
(2010). It has been demonstrated by these and several other authors that, a 
shock-confined shell is unstable to instabilities arising out of shock-induced 
turbulence within layers of the shell. Turbulence leads to enhanced transfer 
of momentum in different regions of the shell that makes it unstable to the 
so called thin shell instability (TSI) (Vishniac 1983, 1994), and raises the 
effective local sound-speed. The stability of the shell depends on the critical 
interplay between the gravitational instability (GI) and the TSI. 

A full fledged analytic treatment of the stability of a shock-confined slab 
is rather complex as has been demonstrated by Vishniac (1994), for instance, 
who showed that such slabs were likely to be unstable to the so called non- 
linear thin shell instability (NTSI). Traditionally, stability analysis of shocked 
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shells and/or slabs, for all practical purposes, have been simplified by excising 
shock dynamics in favour of a simple high-pressure approximation thereby 
eliminating perturbative effects of the TSI, and its non-linear excursion (e.g. 
Whitworth et al. 1994b). Anathpindika (2009, 2010) numerically showed 
that a shocked slab and/or shell is unstable to the TSI, which, soon af- 
ter its formation, develops wiggles on its surface and grows non-linearly. 
Ehlerova & Palous (2002) derived the critical density for a shell to become 
gravitationally unstable. Under a sinusoidal approximation for these per- 
turbations, Anathpindika (2010) deduced an expression for the wavenumber 
of the fastest growing mode. Below, we propose to test the validity of this 
expression for HII shells driven by typical candidate stars listed in Table 5.3 
of Spitzer (1978), and deduce a mass function for clumps condensing out of 
these shells. The case of a shell driven by a supernova blast wave will be 
considered separately. We shall demonstrate that the mass function so de- 
rived is consistent with that reported by Fukui et al. (1999), Yamaguchi et al. 
(2001), and Roslowsky (2005) for massive clouds. In §2 we shall deduce our 
set of equations and demonstrative calculations, including the calculation of 
a mass spectrum for fragments, will be undertaken in §3. We conclude in §4. 



2. Semi-analytic deduction 

2.1. HII shells 

Let us consider a typical source of ionising radiation that emits Nl v c 
number of photons per second. If n p , n e are the respective number of protons 
and e~ per unit volume, then xn e n p a^ is the number of electrons captured 
per cm -3 in the ground state, and the flux of photons flowing through a shell 
of radius, r s , is- 

47T 

—r 3 s xn e n p a {2) = N LyC , (1) 

where is the recombination coefficient that excludes electron captures to 
the ground state, and defined as 

m 2.06 x 10- n Z 2 , , , _ 3 _! 
a = fTT2 <M/3)cm 3 s \ 

at a temperature T; Z is the ionic charge and 02 (/9) is the recombination 
coefficient function corresponding to (see Table 5.2, Spitzer 1978). The 
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ionisation fraction of hydrogen, x, has been set equal to unity so that the 
Stroemgren radius is defined as - 

3 = 3N LyC m 

where the condition of approximate charge neutrality forces n e ~ n p = n. 

The temporal evolution of the ionised shell can be obtained by rewriting 
Eqn. (2) as 

dt Awrfn Airnrf LyCK ' 3 1 J ' 
integration of which yields, 

rj = r s 3 [l - exp(-W 2 ^)]. (3) 

Maximising this equation gives the timescale, transform, over which the shell 
of ionising radiation makes a transition from the initial rarefied phase, to the 
dense phase, which is 

^transform ^ ^y- (4) 

The radius of the expanding shell in the dense phase is given by 

<M = ( l + l«Jii£f, (5) 

r s V 4 r s J 

(e.g. Shore 2007). For an approximately spherical shell, its average volume 
density, p s , is 

_ 3M s hell _ 3M s hell ,„s 

Ps ~ ~4^f~ ~ 4irr 3 s (l + 7a HII t/4r s y 2 / 7 ' 1 ) 

After a time t = transform, the shell cools down to an equilibrium tem- 
perature, T eg , defined by Eqn. (8) below. To estimate the temperature, T e? , 
we shall first account for the likely heating and cooling mechanisms. A cru- 
cial contributor towards heating the shell is the photoionisation of H2, the 
corresponding rate of heating is 

nri~n 2 0(TdkBTi ergs cm" 3 s" 1 ; (7) 
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(Tielens 2005), where fifc = 10 4 K) ~ 1.6 x 10~ 13 cmV 1 , is the recombi- 
nation cooling coefficient. The cooling rate due to collisional excitation of 
carbon is 

n 2 A c ~ (3 x 10' 2 > 2 ( L4 ^ Q _ 4 ) exp (^) ergs cm" 3 ^ 

(Tielens 2005), where Ac ~ 2.6 x 10~ 4 , is the abundance of carbon in the 
ISM. The equilibrium temperature of the gas, T eq , within the shell can be 
estimated using the condition for thermal equilibrium, nTi = n 2 A c , and 
excess energy will be radiated away so that (nTi — n 2 A c ) ~ 3/csT eg /2, i.e. 

2(nr t -n 2 A c ) 

w B — ' (8) 

which is roughly 10 K. 

The temperature, Tj, of the ionised gas is Tj ~ M 2 ^^ Tism, where 
Tism = 100 K, is the average temperature of the preshock ISM and M. is the 
Mach number of the propagating shock, M 2 ~ (vi/a) 2 , a 2 = kj3 ^ M ■ The 
density of the shocked gas, P2 ~ -A4 2 Pi ~ M 2 (n/;%), rhn is the average 
mass of atomic hydrogen. Anathpindika (2010), by performing a perturbative 
analysis, derived an expression for the fastest growing unstable mode in a 
shocked shell, and the wavenumber, k, of this mode is 

k = 7 f E ' 0) 

a 2 _ V2 M _ r °ut 1 
/92 \ r ' m J . 

Here r OMt and r in are respectively the outer, and the inner radii of the shell, 
and P2/P2 ~ (dr-i/dt) 2 = \/ 2 . From Eqn. (3) it follows that, 



dr-j r 3 na' 2 ' 
rft 3 r 



'£^_exp(-W 2 )t). (10) 



For a shell of thickness, dR s , r out = ri n + dR s , so that 1 — ~ dRs , and 

Tin ^"in 

since M sfee// ~ Anr 2 out dR s p s ~ Airr 0Ut dR s Y; s ; S s ~ r ou tp s . For the wavenum- 
ber, /c, calculated using Eqn. (9) above, the corresponding wavelength, 
A c /ump = 27r/A; so that mass of a clump, M c i ump ~ A 2 iump S s , and the number 
of fragments N /rag ~ M sheU /M dump . 

2 f ^ \ ^ 



5 




Figure 1: Velocity characteristics for 5 stars belonging to spectral types between 05 and 
Bl have been plotted. The rightmost characteristic is for a typical 05 star (T^47000 K & 
Nlj,c~ lO 4 ^ 1 ), while the leftmost is for a Bl star (T-22600 K & N LyC ~ lO 4 ^- 1 ). The 
transformation of the shell from R-phase to the D-phase occurs at t ~ transform (Eqn. 
4); thereafter the dense shell progressively decelerates. The velocity of the shell, V s , in the 
R-phase, not plotted here, is at least two orders of magnitude higher. 



and since V s » a , 



4 + v? 



dR, 



1 



1;1 dR x \i 

clump 



M r i umr) ~ V„ 



G 2 (r out p s )' 

The quantity f^^\ ~ 10~ 2 , as will be shown in §3 below, so that 

M dump ~ 10'(^-) 4 ( Pi) ( 10 " 2 ' gCm "> e . (11) 
Vkm/s/ Vr out / V p s / 

Calculations for typical 0-B stars are deferred for §3 below. 
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Figure 2: Plots showing the mass swept up by expanding shells for individual stars. The 
most powerful, 05 star, sweeps up the most massive shell as indicated by the topmost 
characteristic while the characteristics for stars of other spectral type stack below it. The 
temperature of driving stars decreases progressively for plots from top to bottom. 
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Figure 3: The characteristics plotted in green are for the shells driven by hotter O-type 
stars. They have similar thickness-to-radius ratio, the curves to the right, extending up to 
>0.2 Myrs; while the B-type stars which are much cooler drive comparatively smaller shells 
and so the ratio is a little higher, at least by 30%, curves to the left of the former type. 
The integration for these plots was terminated when respective shells became subsonic. 
The corresponding ratio for a SNR is much smaller, implying considerably thinner shells. 
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2.2. Supernova driven shells 

We now discuss the stability of a shell driven by blast waves originating 
from a supernova that injects energy, E s , into the ISM. For a shell having 
average density p s , and radius, R s , 

e.~ ■ (12) 



3 srb \ dt 2 

Integrating Eqn. (12) yields the well-known Sedov solution for the shell 
radius, 

R s {t)=C'{E s /p s ) 1 'H 2 l\ (13) 

where C is a numerical constant of order unity 0. The post-shock tempera- 
ture, Tj, can be calculated using the usual Hugoniot shock conditions for an 
ISM at a preshock temperature, Tism = 100 K, as before. The pressure due 
to the blast-wave, pb ~ pik^Ti/rhn, under the assumption of approximate 
isothermality; all symbols have their usual meaning. The average density 
of the shell, p s ~ A4 2 pi ~ .M 2 (n#m#). Using Eqn. (13) above, we can 
estimate the timescale, ti SOt hermai, over which the shell is likely to acquire 
its equilibrium temperature, T eq , defined by Eqn. (8) above, whence it may 
also slow down substantially, to a sonic or possibly, even sub-sonic speed, 
a = y/kBTisM/mH, is the local sound-speed. 

Then at t tisothermali 

dR s = (^s\ 1/5 .-3/5 = 5a 
dt \p s J isoth ~ l 2 



=> t isothermal ~ 0.217a- 5 / 3 (^) 1/3 . (14) 

VPs / 

The surface density of the shell, S s , is calculated as before and the mass of 
the shell, M sheU 

M shell ~ 10 2 ( n.) 2 ( ( P' ) M ; (15) 



vpc/ VO.lpc/ VlO 21 gcm 3 . 

r\ = r s (t = t isotherma i) , is the radius of the shell at that epoch; this expression 
is also used to calculate the mass of the HII shell. The fastest growing 
unstable mode, as before, is calculated as before, using Eqn. (9) above. The 
calculations for a typical supernova remnant (SNR) are demonstrated in the 
following section. 



L C - (15/Att) 1 / 5 
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3. Results and discussion 

The collect and collapse model : The HII shell 
For illustrative purposes we consider a typical particle density, n ~ 10 cm -3 . 
The minimum timescale for radiative cooling of the shell, t trans f orm , calcu- 
lated using Eqn. (4) is approximately 0.1 Myrs for a typical 05 star; we 
have adopted physical parameters defined in Table 5.3 of Spitzer (1978). 
The shell, at this epoch, enters in to the dense phase whence it acquires 
substantial mass during the snow-plough phase, and undergoes substantial 
deceleration as can be seen from the characteristics plotted in Figs. 1 and 
2. Relatively cooler stars of intermediate mass drive weaker ionisation fronts 
and collect considerably lesser mass. The uppermost characteristic in Fig. 1 
shows that a typical 05 star may drive a shell that has a typical mass of a few 
times 10 6 M . An important parameter associated with the stability of an 
HII shell is its thickness. By a thin shell we imply, j^ 4 - <C 1. It can in fact, 
be easily demonstrated that the ratio of the shell-thickness to radius, may 
rapidly asymptote to a value less than unity. The average surface density of 
the shell, 

V _ r iP g 
— , 

Tout 

which following a little manipulation, leads us to 

dR s Tip s 

r^j _ 

For an HII shell, oc t 4 / 7 , and p s oc t~ 12 / 7 , then for a shell of constant surface 
density, — oc t~ 8 / 7 . For a typical 05 star, this ratio is ~ 0.18. The green 
curves in Fig. 3 shows time evolution of the HII shell thickness-to-radius ratio 
for the set of stars used for demonstrative calculations above. These plots 
suggest, the driven shell is generally thin irrespective of the spectral type of 
driving star. This is consistent with that reported by the GLIMPSE survey 
comprising 600 HII shells (Churchwell et al. 2006). Both, theoretical (e.g 
Vishniac 1983, 1994), and numerical (e.g. Anathpindika 2009, 2010) work has 
demonstrated the susceptibility of thin shells to various shearing instabilities, 
and particularly, to the thin shell instability (TSI). Thin, shock confined 
shells tend to show a greater proclivity towards the TSI that apparently 
dominates the classical Jeans instability. The thermal Jeans mass, M Jeans , 
is 
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which is roughly an order of magnitude smaller than the typical mass of a 
clump, M c i ump , defined by Eqn. (11) above which suggests, perturbations 
triggered by dynamical instabilities in driven shells and/or shocked slabs fa- 
cilitate concentration of material in perturbed regions, the local maximas or 
minimas. This was explicitly demonstrated by Wuensch & Palous (2001). 
The clump on gaining sufficient mass condenses out, which is the essence of 
the collect and collapse model. 

A supernova shell 

For a demonstrative calculation, let us consider a typical supernova that 
releases energy, E s ~ 10 51 ergs, into the ISM that as in the previous case, 
is assumed to have n ~ 10 cm" 3 . The timescale, ti SOt hermai, defined by Eqn. 
(14) above, for a shell of average density ~ 10~ 21 g cm" 3 is ti SOt hermai ~ 77 
Myrs, and the velocity of the shell, V s , 



v, dR 



' s 



dt 



(17) 

*• ''isothermal 



which is ~ 1.5km/s. The corresponding radius of the shell is then, 
R s (t) ~ V s (t)(t = ti. so thermai) ~ 0.11 kpc. The mass of the shell swept up 
immediately follows from Eqn. (15), and M s h e a ~ 10 9 M Q . Similarly, the 
mass of a typical fragment calculated using Eqn. (11) is, M c i ump ~ 10 6 M . 
It might be interesting to calculate the efficiency, 5 with which the initial 
energy, E s , is converted into mechanical energy whence matter in the ISM is 
swept up. This efficiency is defined as, 

5 = 1- ( Qa ^ reMated ) ~ 0.999, (18) 

^ a SNR ' 

which suggests that a large proportion of the initially injected energy is lost 
in heating the ISM. The attenuated sound speed, a attenuated defined by Eqn. 
(19) below accounts for the effects of turbulence generated by various hydro- 
dynamic instabilities which were not directly included in the perturbative 
analysis that led to Eqn. (9), and 



(^attenuated 



P2(l_ T 'out 
P2 



(19) 



The GLIMPSE survey of supernovae shells produced a catalogue of 95 
SNRs (Reach et al. 2006), none of which have shown unambiguous evidence 
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in favour of the collect and collapse model, although a sample of SNRs have 
shown association with OH masers. However, the source of these masers 
is unclear and could perhaps be due to sporadic star-formation triggered in 
molecular clouds over-run by SNRs. Despite this being the case, one may 
assert that possible condensation in SNRs may take a while, as suggested 
by the magnitude of t isot hermah which is roughly two orders of magnitude 
larger than the timescale over which the HII shell produced fragments. The 
resulting fragments are therefore likely to be at least as massive as those 
forming in the former case. 

3.1. The clump mass spectrum 

The number of clumps, N dump) condensing out of a shell is N dump ~ 



This equation defines the number of clumps in an interval (k, k + dk), the in- 
tegration of which, over the wave number domain yields the mass spectrum 
for putative clumps. Figure 4 shows a typical mass spectrum for clumps 
condensing out of HII shells, obtained via a Monte-Carlo integration of Eqn. 
(20). The integration was performed for 10,000 realisations of shell fragmen- 
tation, each producing N dump number of clumps. Thus, we have ~180,000 
clumps for the 05 star and fewer, ~ 150, 000, for an 09 star. The relatively 
smaller shell driven by the latter star not only produces fewer fragments but 
also, those that are comparatively less massive than those in the former case. 
However, we note that extremely large clumps, M c i ump /M Q > 10 6 , may only 
be sparingly produced as is evident from the mass spectrum in Fig. 4. The 
mass spectrum for stars of either spectral type is similar and a power-law, 
dN / dM c i urnp oc M~ t ^ , appears to fit the derived spectrum reasonably well 
in either cases; (3 ~ 1.6, and 1.5 respectively. The spectrum derived here is 
consistent with the one obtained for large clouds by Fukui et al. (1999), and 
Roslowsky (2005). 



M s hell 
M clump 



. Then 




or equivalently, in terms of the wavenumber, k, 




(20) 
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Similarly, we also obtain the mass spectrum for fragments condensing out 
of the SNR driven shells. The spectrum for this case is plotted in Fig. 5, 
and a power-law similar to the one in earlier two cases, fits the spectrum 
reasonably well. In fact, this latter spectrum is similar, j3 — 1.6, to that for 
the fragments resulting from the HII shell driven by an 05 star; spectrum for 
the SNR was derived for 10 5 realisations of shell fragmentation, and shows 
considerable shift towards a higher mass in comparison to the distributions 
shown in Fig. 4. This is probably because the velocity of an HII shell in 
its dense phase decays more rapidly than that of an SNR, which pushes up 
the attenuated sound-speed, a a ttenuated, for the latter that in turn raises the 
mass, M dump , as seen in Fig. 5. 

Star forming clumps as massive as ^2500 M Q have been reported in the 
HII shell N49 and RCW34 (e.g. Zavagno et al. 2010, Bik et al. 2010). A 
similar treatment of the problem was presented by Whitworth et al. (1994a,b) 
who arrived at a minimum clump mass of only a few tens of M Q , which in 
the light of present findings, appears somewhat conservative. Our claim 
more massive clumps is supported by recent observations described above. 
It is well known, the evolution of a cloud is governed by the dynamical 
effects associated with the complex interplay between self-gravity, and other 
contributing factors such as turbulence and the magnetic field that support 
a clump against the former. It is therefore crucial to predict clump masses 
with reasonably good accuracy; under-estimation of the mass of putative 
clumps, for a given radius, which in the present case is roughly equal to the 
wavelength of the unstable mode, will lower the average density and thus 
raise the clump-lifetime. An increase in the longevity of clumps will also 
possibly make events such as clump-clump collisions more probable, than 
they are known to be. Inter-clump collisions, according to simulations, could 
lead to bursts of star-formation with fewer low-mass stars (e.g. Chapman et 
al. 1992, Anathpindika 2009), that will tend to shift the stellar initial mass 
function towards a top-heavy distribution, away from the widely reported 
lognormal form. 

4. Conclusions 

We have examined the stability of thin, dense shells driven by powerful 
ionising radiation originating from massive, early type stars and/or blast 
waves from a supernova. Our work here shows that fragmentation of these 
shells is likely to produce large clumps, with masses typically > 10 3 M . 
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Figure 4: The clump mass spectrum generated via a Monte Carlo integration of Eqn. 
(20). The green spectrum for a typical 05 star peaks at ~ 10 3 M0; the peak shifts a 
little rightward, closer to ~ 10 4 M © for an 09 star that is much cooler (red histogram). 
Formation of massive clumps > 10 6 Mq is possible, though only sparingly. A power-law 
fit of the type dA f N — ocM^ agrees reasonably well with the derived spectrum. Also 
see text below. 
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Figure 5: The clump mass spectrum for fragments condensing out of an SNR. The general 
features of this spectrum are similar to those of the spectrum for clumps produced via 
fragmentation of a HII shell. See text below. 
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This calculated range of fragment masses is consistent with that reported 
via observations of HII shells, and the fragments so formed could spawn a 
second generation of stars, which may in turn trigger the next generation 
of stars in the surrounding ISM. The simple calculations discussed above 
appear to suggest that the ratio of thickness-to-radius for a shell evolves 
only weakly with time, and that calculated here (see Fig. 3; < 10" 1 ), is 
consistent with the values reported in the GLIMPSE survey for HII shells. 
Similar results are obtained for an SNR, which, however, appears to be much 
thinner compared to the shells driven by OB stars. This could perhaps be the 
reason why SNRs often appear filamentary in the IRAC bands (e.g. Reach 
et al 2006). Simulations studying the evolution of thin shells such those by 
Anathpindika (2010) have demonstrated the dominance of TSI in shocked 
shells; the surface of the shell was also shown to develop ripples, similar to 
breathing modes on fluid surfaces (Fig. 3 in Anathpindika (2010)). While 
wiggles on the shell surface associated with the TSI though coplanar, are 
generally orthogonal to its surface, and so in case of a shell in the plane of 
the sky, cannot alone account for the reported filamentary nature. However, 
breathing modes coupled with local magnetic field could perchance explain 
the occurrence of filaments reported by Reach et al (2006). 
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